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One of the most striking universal properties of the high-transition- 
temperature (high-T c ) superconductors is that they are all derived from the 
hole-doping of their insulating antiferromagnetic (AF) parent compounds. From 
the outset, the intimate relationship between magnetism and superconductivity 
in these copper-oxides has intrigued researchers []I|-f§]. Evidence for this link 
comes from neutron scattering experiments that show the unambiguous pres- 
ence of short-range AF correlations (excitations) in cuprate superconductors. 
Even so, the role of such excitations in the pairing mechanism and superconduc- 
tivity is still a subject of controversy fl5]. For YBa 2 Cu 3 06+:r, where x controls the 
hole-doping level, the most prominent feature in the magnetic excitations spec- 
tra is the "resonance" [pHllfl- Here we show that for underdoped YBa 2 Cu 3 6 6 , 
where x and T c are below the optimal values, modest magnetic fields suppress 
the resonance significantly, much more so for fields approximately perpendicular 
rather than parallel to the CuC>2 planes. Our results indicate that the resonance 
measures pairing and phase coherence, suggesting that magnetism plays an im- 
portant role in the superconductivity of cuprates. The persistence of a field 
effect above T c favors mechanisms with preformed pairs in the normal state of 
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underdoped cuprates p2|-|14|. 

Neutrons carry magnetic dipoles which are scattered by the electron spins in the sample 
with a probability, dependent on the neutron direction and energy, directly proportional to 
the Fourier transform in space and time of the two-spin correlation function, S(Q,u). The 
wave-vector Q and energy Uu are simply the difference between the momenta and energies, 
respectively, of the incident and scattered neutrons. YF^CusOe+x [denoted as (123)06+a;] 
with x of 0.6 was selected because it is a thoroughly characterized high-T c superconductor of 
which it is possible to fabricate the large single crystals required for neutron measurements. 
Equally important, thermodynamic data indicate that even modest fields have a large effect 
on the electronic entropy in the superconducting phase [jn^[nj|. Thus, if the electronic 



entropy changes around T c are mostly due to spin excitations |TT|, they must respond to an 
external magnetic field in qualitatively the same way [17] . 



(123)06+a; contains pairs of Cu02 layers separated by charge reservoir layers. We label 
wave vectors Q = (q x ,q y ,Qz) m units of A as (H,K,L) = (q x a/27i,q y b/27i,q z c/27i) in the 
reciprocal lattice units (rlu) appropriate for the orthorhombic unit cell, for which the lattice 
parameters are a = 3.83, b = 3.88 (~ a), and c = 11.743 A if x = 0.6. For x — 0, (123)06+^ 
is an antiferromagnetic (AF) insulator where the spin on each Cu 2+ ion is antiparallel to 
its nearest neighbors both within as well as in the immediately adjacent layer ||18|| . As the 



oxygen content x is increased, the material becomes a metal and eventually, near x — 1, an 
excellent high-temperature (T c ~ 93 K) superconductor. At x near 0.6, (123)06+x is also 
a good superconductor with T c reduced to ~60 K. During the evolution from insulator to 
superconductor, static antiferromagnetism is lost while magnetic (spin) fluctuations centered 
at the x = AF Bragg positions, often referred to as Qo = (n, ir) (refs J^^TTJ] ) , persist. For 
highly doped (123)0 6+:r , the most prominent feature in the spin fluctuation spectrum is 
a sharp resonance which appears below T c at an energy of 41 meV (refs When 
scanned at fixed frequency as a function of wave vector, the sharp peak is centered at 
(11,71) (refs [§]-|8|]) and its intensity is unaffected by a 11.5 T field in the ab-plane [fTS 



In our underdoped (123)06.6 (T c = 62.7 K) (ref. ||), the resonance occurs at 34 meV 



and is superposed on a continuum which is gapped at low energies 111]. For frequencies 



below the resonance, the continuum is actually strongest at a quartet of incommensurate 



positions disposed symmetrically about (71*, 7r) (refs pO|j2T|| ). Both the resonance and the 
incommensurate fluctuations appear in the acoustic channel for the bilayers, meaning that 
they are due to spin correlations which are AF between the planes. The corresponding 
neutron scattering intensities exhibit a modulation sin 2 (7r<iL/c), where d (3.342 A) is the 
spacing between the nearest-neighbor CuC>2 planes, along c. In contrast, optical fluctuations 
are associated with ferromagnetic correlations between planes and vary as cos 2 (irdL/c) along 
c. The optical fluctuations for x = 0.6 have a gap of ~50 meV, well above the acoustic gap 
of 20 meV (ref. 0). 

We mounted the crystal in three different orientations inside a split coil (to provide 
neutron beam access) 7 T vertical field superconducting magnet to establish the effects of 
a field approximately along the c-axis (B ~ ||c-axis, Fig. la) on the acoustic, along the 
c-axis (5 1 1 c-axis) on the optical, and in the ab-plane (5||a6-plane, Fig. Id) on the acoustic 
fluctuations. Raw neutron scattering data for (123)0 6+x in the large frequency range covered 
here contain considerable non-magnetic background. The background, however, is weakly 
temperature-dependent for T < 70 K, which makes temperature difference spectra a useful 
guide to the magnetic contributions. Figure lb shows the effect of superconductivity on 
S(Q,lj) at Q = (1/2,3/2,1.8) in the form of difference spectra between 10 K (< T c ) and 
70 K (T c + 7.3 K). At zero field (B = T), we reproduce previous results in both the 
B ~ I |c-axis (Fig. lb) and the B\ |afe-plane (Fig. le) sample orientations. The data reveal the 
positive peak at 34 meV due to the increased resonance amplitude below T c . Imposition of a 
6.8 T field nearly along the c-axis results in the qualitatively similar temperature difference 
spectrum (Fig. lc). However, the intensity gain of the resonance is suppressed by ~30% 
compared to that at zero field (Fig. lb). This is remarkable because the magnitude (6.8 T) 
of the applied field is much less than the upper critical field [B c i = 45 T) for the material 



Tfifl , and the corresponding Zeeman (magnetic) energy is, at igfiB (~ 0.8 meV assuming 



g = 2 and S = 1), much smaller than the 34 meV resonance energy and the thermal energy 



ksT c (~6 meV). 

In common with other layered superconductors, the thermodynamic properties of the 
cuprates are strongly anisotropic with respect to the direction of the applied magnetic field 
|T6| . A field perpendicular to the conducting planes generally causes a stronger suppression 
of the superconductivity than one applied within the plane. To see whether the resonance 
intensity depends similarly on field direction, we imposed a 6.8 T field parallel to the ab- 
plane (Fig. If). A noticeable, but much smaller suppression (< 10%) of the resonance is 
observed. Thus, the anisotropy of the field effect on the resonance reflects the anisotropy 
of other physical superconducting properties. This is the clearest evidence to date that the 
resonance is very sensitive to the superconducting pairing. 

Although the constant-Q difference spectra (Fig. 1) are excellent for determining the 
intensity gain of the resonance from the normal to the superconducting state, such measure- 
ments do not provide information about how the field affects the momentum distribution of 
the excitations. We therefore scanned wave vector in the B ~ ||c-axis geometry (Fig. la) at 
energy transfers below, at, and above the resonance. At 150 K, the constant-energy scans 
at the resonance energy (Fig. 2a) show a broad peak centered at (it, it) on a sloped linear 
background. There are no observable intensity differences in the profiles for B = and 6.8 
T. On cooling to 70 K (Fig. 2b), the resonance narrows in width and grows in intensity, 
but is also clearly suppressed at 6.8 T. At 10 K (Fig. 2c), the suppression of the resonance 
at 6.8 T is more severe. Direct inspection and fits of Gaussian profiles to the data (Fig. 
2, a-c) indicate no appreciable field-induced change in the momentum dependence of the 
resonance peak. In contrast, the thermal fluctuations associated with warming from 10 K 
to 70 K yield a decrease in Q-integrated intensity which is numerically similar to that found 
on application of a 6.8 T field at 10 K, and makes the profile considerably broader. 

Below the 20 meV spin gap in the superconducting state [|1TJ, the constant-energy scan 
at 16 meV in the 6.8 T (Fig. 2d) is featureless, that is, there are no strong field-induced 
AF correlations at this energy. At 24 meV, where incommensurate spin fluctuations were 
observed |20|j2l[| the field seems to suppress the constant-energy profiles in both the normal 



(Fig. 2e) and the superconducting state (Fig. 2f). The effect is more difficult to see 
than at the resonance energy although it could represent a similar fractional change in a 
considerably smaller signal above background. Finally, when we move to an energy just 
above the resonance, we discern no field-induced change in the Q-dependent profiles at both 
70 K and 10 K (Fig. 2g and h). 

To demonstrate the effect of a magnetic field on the energy dependence of the excitations 
in (123)0 6 . 6 , we show the experimentally measured difference spectra, S(Q ,uj,B = 0T) — 
S(Q ,lu,B = 6.8T), for B ~ ||c-axis (Fig. 3, a-c) and in the ab-plane (Fig. 3d and 
e) at different temperatures. Since one expects the lattice contribution (phonons) to the 
scattering to be field-insensitive, this procedure will yield the net change of the magnetic 
intensity induced by the field. At 150 K, a 6.8 T field has no effect on the spin fluctuations 
in the measured energy range for B ~ ||c-axis (Fig. 3a). On cooling to 70 K, the difference 
spectra show a broad peak for B ~ ||c-axis (Fig. 3b) and no discernible feature for B\\ab- 
plane (Fig. 3d). Higher resolution measurements in the B ~ ||c-axis geometry at 6.2 T field 
(Fig. 3f) confirm the result of Fig. 3b. In the superconducting state (Fig. 3c, e, and g), the 
difference spectra show a sharp peak at 34 meV for both field orientations, again confirming 
that the resonance is the feature most obviously affected by the magnetic field for spin 
fluctuations below 50 meV. In addition, the sharpness here contrasts with the diffuse nature 
of the difference just above T c , thus showing the sensitivity of the resonance width, or decay 
rate, to true bulk superconductivity even in a material with an extended pseudogap regime 
above T c . Simple summation as well as least-squares fits to the data (Fig. 3f and g) indicate 
that the energy-integrated effect of the magnetic field just above T c is comparable to that 
well below T c . The narrowing of the widths in the field-induced difference spectra together 
with the approximate preservation of weight in these spectra suggest that the resonance 
width depends strongly on the phase coherence time, while its integrated weight measures 
the pair density. 

In Fig. 4 we plot the temperature dependence of the neutron scattering intensity at the 
resonance energy at zero and 6.8 T field. Most growth in the field-dependent peak signal 
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tracks the growth in the resonance itself on entering the superconducting state. The inset 
shows the field dependence of the normalized intensity. The solid line is a linear fit to the 
data, i.e., //% = 1 — B/B c h ar with B c h ar = 36 T. Because B c h ar is not far from the upper 
critical field B& = 45 T estimated for our sample, it is reasonable to associate the sharp 
resonance peak with the superconducting volume fraction which for simple superconductors 
is proportional to 1 — B/B c2 . 

Figures 2 and 3 indicate that wherever we look at fixed temperature, the external field 
either has no effect or reduces the magnetic scattering signal. The total moment sum rule 



22| , |23|| states that the magnetic structure factor, when integrated over all wave vectors and 
energies should be a temperature- and field-independent constant. Therefore, compensatory- 
increases in signal should occur elsewhere in Q-u space. Possible destinations for the lost 
spectral weight include elastic scattering, new or softened optic modes, a subgap contin- 
uum, and other continuum scattering spread sufficiently over Q and u space so as to be 
unobservable at any particular Q and to. Searches for the most obvious elastic scattering, 
namely that due to the proposed field-induced AF order in cuprate superconductors pi 



have been fruitless. Scans in the [H, K, 0] plane with B\ |c-axis revealed no evidence for optic 
mode scattering at hu = 34, 41, and 50 meV apart from that already present at zero field. 
In addition, we have not yet found field-induced subgap scattering (Fig. 2d). While our 
searches for the missing spectral weight were by no means exhaustive, a 6.8 T field does not 
induce obvious new excitations around (jr, it) with energies between 10 and 50 meV that are 
big enough to compensate the spectral weight loss of the resonance. It is noteworthy that 
similar sum rule violations also exist for spectra collected using angle resolved photoemission 
|25| , scanning tunneling microscopy |26) , and infrared spectroscopy |27| in some high-T c su- 



perconductors. However, we caution that an explanation of the superficially similar effects 
in these measurements must take into consideration that different experimental technique 
probe quite different correlations. 

We have discovered that a very modest field applied to (123)06.6 yields a very significant 
reduction in the 34 meV resonance. Therefore, the resonance can be reduced either by 
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warming or applying a magnetic field, with the highest amplitude and narrowest profiles 
both in energy and wave vector occurring in the bulk superconducting state. Even more 
surprising is the persistence of the effect above T c (Figs. 2-4). Although there are no explicit 
microscopic predictions about the effects of magnetic fields on the resonance, our finding has 
important implications for understanding the mechanism of high-T c superconductivity. In 
particular, the larger c-axis field effect presents another challenge to the interlayer tunneling 
theory, according to which the resonance should be strongly affected by the in-plane field 
that disrupts the coherent Josephson coupling along c-axis p8 |. On the other hand, the 
data are consistent with mechanisms where the dominant loss of entropy on entering the 
superconducting state is due to the growth of magnetic correlations in the Cu02 planes. A 
modest field reduces the resonance just as it reduces the specific heat anomaly |1| near T c , so 
that the temperature derivative of the resonance intensity tracks the specific heat not only as 



the hole density is varied |TTJ, but also as the magnetic field is imposed for fixed x (ref. |fL7|j ). 
The anisotropic field effect reflects the anisotropy of other physical properties such as the 
magnetic penetration depth and coherence length. Thus, the persistence of the effect above 
T c would naturally be due to superconducting fluctuations or incoherent pairing which also 
display this anisotropy in the normal state. Such fluctuations exist in models postulating a 
Bose-Einstein condensation of pre-formed pairs at Tc in the underdoped cuprates fl2T-|i~4|l . 
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FIGURES 

FIG. 1. Reciprocal lattice diagram and neutron scattering results. For the experiment, we 
used the HB-3 triple-axis spectrometer at the High-Flux Isotope Reactor of Oak Ridge National 
Laboratory. The collimations were, proceeding from the reactor to the detector, 48-40-60-120 
min (full-width at half-maximum) and the final neutron energy was fixed at either Ef = 30.5 
meV or 14.78 meV. The monochromator, analyzer and filters were all pyrolytic graphite. The 
low temperature data were collected after field cooling the sample from 200 K. a, To determine 
the effect of a field along the c-axis on the acoustic spin fluctuations, we first align the crystal in 
the [H, K, 0] zone. The sample is then rotated a degrees around the [3H, —H, 0] direction, where 
a = arctan[(2vrL/c)/(27r v /^ 2 + (3H) 2 /a)] « 20.6° for H = 0.5 and L = 1.8 rlu. such that the 
horizontal scattering plane is spanned by wave vectors along the [1, 3, tan aVWc/a] and [3, —1,0] 
directions. In this geometry, we can perform both constant-Q (Fig. lb and c) and constant-energy 
(Figs. 2, a-h) scans at the acoustic mode intensity maximum (L = 1.8 rlu) while still keeping 94% 
of the applied (vertical) field along the c-axis (B ~ || c-axis). To see the influence of a field along 
c on the intensity of the optical mode, we oriented the c-axis of the crystal perpendicular to the 
scattering plane (B\ | c-axis). This geometry allows any wave vector of form (H, K, 0), that is, at 
the L = intensity maximum of the optical mode, to be reached, d, Conventional [H, 3H, L] zone 
geometry where the applied field is in the ab-plane along the [3,-1,0] direction, b, c, Difference 
spectra of the neutron intensities between T = 10 K (< T c ) and T = 70 K (T c + 7.3 K) at wave 
vector Q = (0.5, 1.5, 1.8) for zero and 6.8 T field in the B ~ || c-axis geometry, e, f, Similar data in 
the .B 1 1 aft-plane geometry. Because the differences in sample geometries, the actual intensity gain 
of the resonance in the aft-plane experiment is slightly higher than that of the B ~ || c-axis case 
at zero field for the same incident beam monitor (300 counts which correspond to ~8 min per point 
at hu = 34 meV). We normalized the intensity gain of the resonance in these two experiments at 
zero field and multiplied the same scale factor to the finite field case. 
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FIG. 2. Effect of a magnetic field on momentum dependence of the spin excitations. Con- 
stant-energy scans of the neutron scattering intensity as a function of wave vector in the B ~ | |c-axis 
geometry, a-c show the results at the excitation energy of 34 meV; d, 16 meV; e, f, 24 meV; and 
g, h, 41 meV. The open and filled circles represent data for identical scans at zero and 6.8 T field, 
respectively. The dashed and solid lines are Gaussian fits to the data. 

FIG. 3. Effect of a magnetic field on the energy dependence of the spin excitations. Difference 
spectra of the neutron scattering intensities between zero and finite magnetic field in a-c and d-e 
show the low-resolution results obtained with Ef = 30.5 meV at wave vector Q = (0.5,1.5,1.8) 
in the B ~ ||c-axis and B||a6-plane geometries, respectively. High-resolution measurements in the 
B ~ ||c-axis geometry with Ef = 14.78 meV are shown in f-g. The scattering should be centered 
around zero as shown by the solid lines in a and d if spin fluctuations are not affected by the 
field. The solid lines in c and e-g are unconstrained least-squares fits to the data by Gaussians 
on floating constant backgrounds. The solid line in b is a guide to the eye. Positive scattering 
(above zero) indicates a net loss in the acoustic magnetic spectral weight when the field is applied. 
Strong phonon scattering at ~20 meV in the raw constant-Q scans yielded the large error bars in 
the difference spectra around that energy. 

FIG. 4. Effect of a magnetic field on the temperature dependence of the resonance and the 
field dependence of the resonance at ~10 K. The data are collected in the B ~ ||c-axis geometry. 
The open and filled circles show results at Q = (0.5, 1.5, 1.8) and Huj = 34 meV for zero and 6.8 T 
field, respectively. The pseudogap temperature T* and T c of the sample are indicated by arrows. 
The solid and dashed lines are guides to the eye. The inset shows the field dependence of the 
normalized resonance intensity at approximately 10 K. The filled squares are high-resolution data 
at Q = (0.5, 1.5, 2) with E f = 14.78 meV. The solid line corresponds to I/I = 1 - (B/B char ) with 
Bchar = 36 T while the dashed line represents I/Iq = 1 — (B / B c h a r) 1 ^ 2 , where B c h ar = 208 T. 
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